arXiv: 1505.0667Ivl [math.DG] 25 May 2015 


Singularities of the geodesic flow on surfaces with 
pseudo-Riemannian metrics 

A.O. Remizov* and F. Tari^ 

May 26, 2015 


Abstract 

We consider a pseudo-Riemannian metric that changes signature along a 
smooth curve on a surface, called the discriminant curve. The discriminant 
curve separates the surface locally into a Riemannian and a Lorentzian domain. 

We study the local behaviour and properties of geodesics at a point on the 
discriminant where the isotropic direction is tangent to the discriminant curve. 

II 

1 Introduction 

The work in this paper is a part of an ongoing research on understanding the geometry 
of surfaces endowed with a signature varying metric (see, for example, 0 [ini Fa [H 
[T51 FSl FZl Fa Fa Ea Ell Ea EH Ea Ea Ea ) • We consider here the behavior of geodesics 
on a surface S endowed with a pseudo-Riemannian metric given in local coordinates 
by 

ds'^ = a{x, y) dx^ + 2b{x, y) dxdy + c{x, y) dy'^ (1.1) 

where the coefficients a, b, c are smooth (that means C°° unless stated otherwise) 
functions on an open set t/ C 

We assume in all the paper that the discriminant function A{x,y) = {ac — b‘^){x,y) 
vanishes on a regular curve which is called signature changing curve or simply 
discriminant curve of the metric fll.ip . The discriminant curve St separates (at least, 
locally) the surface S into a Riemannian (A > 0) and a Lorenzian (A < 0) domain. 
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At any point in the Lorentzian domain, there are two lightlike or isotropic directions 
that consist of vectors with zero length and there is one double isotropic direction at 
any point on The isotropic curves are integral curves of the equation 

a{x, y) dx^ + 26(x, y) dxdy + c{x, y) dy"^ = 0. (1.2) 

Isotropic curves are geodesics (except in the case below) in the metric fll.ip when 
dehned, for instance, as extremals of the action functional (the arc-length parametriza- 
tion is not dehned for these curves), see A non-isotropic geodesic is called timelike 
{spacelike) if > 0 (ds^ < 0) along the geodesic. The exception is when the isotropic 
curve coincides with the discriminant curve and is a singular solution of equation 
fll.2p . i.e., it is the envelop of one-parameter family of isotropic curves. This is the case 
Z (Table [T]) and the reason why the singular solution is not a geodesic is explained in 
Appendix B. 

When the unique isotropic direction is transverse to the isotropic curves form a 
family of cusps (conhguration C in Figured]). At points where the isotropic direction is 
tangent to the isotropic curves have generically one of the conhgurations Dg, Dn, Dj 
in Figure HI Finally, the conhguration Z in Figure HI occurs in the case when the 
isotropic direction is tangent to S> at all points on S>. (See Proposition 12.21 for more 
details.) 







.Vlli 

z 


Figure 1: Top hgures represent the integral curves of the held fl2.2p on the isotropic 
surface The bottom hgures represent their projections to S, i.e., isotropic geodesics. 
The dashed lines represent the criminant (in the top hgures) and the discriminant curve 
(in the bottom hgures). 

The geodesic how generated by the metric fll.ip has singularities at every point 
q & St. The geodesics cannot pass through g G ^ in an arbitrary tangent direction 
but only in certain directions said to be admissible, and the isotropic direction is always 
one of them. Singularities of geodesic how at generic points q E S (excluding isolated 
points) are studied in [101 dB [22] • The excluded isolated points are the following: 
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(a) points where the isotropic direction is tangent to 

(b) points of intersection of ^ with the closnre of the parabolic set of S. 

In this paper, we complete the stndy in [101 EU [22] for the case (a). We give in 
Theorems 13.21 and 13.31 the conhgurations of the geodesic in case (a) that correspond 
to the conhgnrations Dg and Dn,Df of the isotropic geodesics, respectively. We start 
onr stndy with the conhgnration Z fTheorem 13.ip . which is not generic. However, it 
provides a graphical illustration of some properties common with the cases Dg, Dn, Df 
(for instance, there are no geodesics outgoing from the point of interest to the Lorenzian 
domain). 

In Appendix A we recall some results on local normal forms of resonant and non¬ 
resonant vector helds that we use in the paper. Finally, in Appendix B we briefly 
consider naturally parametrized geodesics dehned as extremals of the action functional. 


2 Isotropic curves and geodesics 

As the discriminant curve & is supposed to be a regular curve, the coefficients a, b, c 
of the metric fll.ip cannot vanish simultaneously at any point in a neighborhood of S'. 
We assume, without loss of generality, that near the point of interest the coefficient 
c{x, y) > 0 and write equation fll.2p as the implicit differential equation 

F(x, y,p) := c{x, y)p^ + 2b{x, y)p + a{x, y) = 0, (2.1) 

where p = dy/dx is the non-homogeneous coordinate in the projective tangent bundle 
PTS. At any point q = {x,y) in the Lorenzian domain (A < 0), equation fl2.ip has 
two simple real roots which correspond to the pair of transversal isotropic directions at 
q. At any q E S, equation fl2.ip has a double root po{q) = --{q), which corresponds 
to the unique isotropic direction at q. Generically, the direction pq is transverse to S 
at almost all points while tangency (of the hrst order) can occur at isolated points on 
S only. 

In the (x, |/,p)-space, equation fl2.ip dehnes a surface ^ called the isotropic surface 
of the metric fll.ll) . On ^ is dehned a direction held which is the intersection of the 
contact planes dy = pdx with the tangent planes of This direction held is parallel 
to the vector held 

-Fpd^ + -pFpdy - -(F^ + pFy)dp. (2.2) 

Formula fl2.2p dihers from what is generally used by the factor |. This is for conve¬ 
nience in the computations carried out in this paper. 

Although the vector held 02.21) is dehned in the (x, |/,p)-space, we will consider it 
only on the isotropic surface which is an invariant surface of 02.2p . The isotropic 
curves are the projections of the integral curves of 02.2p on ^ to the (x, |/)-plane 
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along the p-direction, called vertical. The set of points on the surface where the 
projection tt: ^ —)■ is singular (i.e., where F = Tj, = 0 ) is called the criminant 
of equation fl2.ip . That is, the criminant consists of the points {q,Po{q)) with q E 
The projection of the criminant on the (x, |/)-plane is the discriminant curve 3i. 

Remark 2.1 In a neighborhood of any point g* G there exist local coordinates 
where the metric is diagonal (i.e., 6 = 0 ) and the discriminant curve ^ is given by the 
equation a(x, y) = 0. Then po(<?) = 0 at all g G the derivatives a^,, ay do not vanish 
simultaneously (since A is a regular function) and c(g*) ^ 0. We assume, without loss 
of generality, that c(g*) > 0 . 

As in Riemannian geometry, non-parametrized geodesics in pseudo-Riemannian 
metrics can be dehned as extremals of the length functional or as extremals of the 
action functional after ignoring the natural parametrization, see [22] . In the hst case, a 
special problem arises. The Lagrangian of the length functional has the form L = \/F, 
where F dehned in fl2.ip vanishes on the isotropic surface and the corresponding 
Euler-Lagrange equation is not dehned on Also, the arc-length parametrization is 
not dehned for isotropic curves. However, such a problem does not arise if we dehne 
geodesics as extremals of the action functional. 

The standard projectivization TS* —?■ PTS, which means forgetting the natural 
parametrization of geodesics) transforms the held fl2.2p into a direction held on PTS 
parallel to 

2Ad^ + 2pAdy + Mdp, (2.3) 

3 

where M is a cubic polynomial in p given by M{q,p) = ^ Pi{q)p^ with the coefficients 

i=0 

Po = a{ay - 2b^) + a^b, 

Pi b(3ay 2b^ axC 2aCx^ 

P2 = b{2by — 3cx) + 2ayC — acy, 

/X 3 c(2by C 3 ;) bCy. 

See |22| for more details. By Theorem 1 in [TO], the isotropic surface ^ is an invariant 
surface of the held fl2.3p . and integral curves of the 2-distribution dy = pdx lying on ^ 
are integral curves of 02.31) . Hence all isotropic curves (except for those that coincide 
with with an appropriate choice of parametrization are extremals of the action 
functional. Consequently, they are geodesics. 

Consider the held 02.31) at a point (g,p), q E S'. Since A(g) = 0, the vertical line 
passing through g is an integral curve of 02.3p . If M{q,p) 7 ^ 0, it is the unique integral 
curve of 02.3p passing through (g,p). The projection of this curve to the (x, |/)-plane 
along the p-direction is the point g, so is not a geodesic. Hence geodesics can pass 


4 










through a point g G ^ only at admissible tangential directions p given by M{q,p) = 0. 
One can show (see M) that 

M{q,p) = ^{p-po){2{A^+pAy) + Mp), VgG^. (2.4) 

From fl2.4p . it follows that the isotropic direction po is always a root of the cubic 
polynomial M{q,p). Furthermore, if p is a double root of M{q,p), then the direction 
p is tangent to the curve ^ at q. Indeed, if p 7 ^ Po, then substituting M = Mp = 0 in 
fl2.4p . we get A^. +pAy = 0. If p = po, then differentiation fl2.4|] with respect to p and 
substitution M = Mp = 0 yields A^ + PoAy = 0. 

Let K denote the Gaussian curvature of (Generically K{q) 00 as q ^ S>.) 

The function Ki = A^K is well defined and is smooth on the whole surface S and is 
an isometric invariant. The set iFi = 0 is the closure of the parabolic set of S. 

Before giving some properties of M, we introduce the following notation: po rfi ^ 
means that the direction po is transverse to the curve S', Poll-^ and ord(po||^) = 1 
mean respectively that po is tangent to 3! and has first order tangency with 3. Finally, 
(£ 1 , £ 2 ) and (Ai, A 2 , 0) denote the spectrums of the linear parts of the vector fields fl 2 . 2 p 
and fl2.3p at the point (g,Po), q & 3, respectively. 

Proposition 2.1 Let q ^ 3, then the following statements hold. 

(i) The cubic polynomial M{q,p) has a unique real root po if and only if Ki{q) < 0. 

(ii) Po is a simple root of M{q,p) if and only if po &\ 3 at q, or equivalently, the 
p-component of the field 02.211 does not vanish at (g,po). 

(hi) M{q,p) has three distinct real roots po,pi,p 2 if and only if Ki{q) > 0 andpo rfi 3 
at q. 

(iv) M{q,p) has a double root Pi = P 2 7 ^ Po if and only if Ki{q) = 0 and po &\ 3 at 
q. The double non-isotropic direction pi is tangent to 3 at q. 

(v) M{q,p) has a double root po = Pi 7 ^ P 2 if and only if Ki{q) > 0 and pq\\ 3 at q. 

(vi) The isotropic direction Pq\\ 3 at q if and only if 02.2p vanishes at (g,po). Then 

+ ^2 7 ^ 0. Moreover, ord (po||^) = 1 if and only if Si, 7 ^ 0, that is, (g,Po) is 
a node, saddle or focus of 02 . 2 p . 

(vii) If Poll 3 at all points on 3, then the p-component of the field 02 . 2 p vanishes on 
the criminant. Hence the criminant consists of singular points of 02.21) with one 
of the eigenvalues £ 1,62 equal to zero and the other is distinct from zero. 

(viii) If {q,p) is a singular point of the field 02.31) . then Ai = 2{A„,{q) -\-pAy{q)) and 
X 2 = Mp{q,p). 
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Proof It is convenient to take the point q to be the origin and choose the local 
coordinates in Remark 12. II Then fl2.4p becomes 

M{0,p) = -cpQ{p), Q{p) = Cxp^ - 2ayp - a^, 

the qnadric polynomial Q has discriminant 6 = + a^. Using the Brioschi formnla 

PI. we get Ad(0) = c(0)(5(0)/4. Since the coefficient c is positive, it follows that A"i(0) 
and (5(0) have the same sign. 

For (i), the cnbic polynomial M(0, p) has only one real root po if and only if (5(0) <0 
or (5(0) = 0 and po = 0 is a donble root of Q{p). In the second case we have 03 ,( 0 ) = 0, 
which together with (5(0) = 0 gives 03 ,( 0 ) = 0 ^,( 0 ) = 0. This contradicts the assnmption 
that A = oc is regnlar at the origin. 

The proof of (ii) follows by direct calcnlation. For (iii), M{0,p) has three distinct 
real roots if and only if Q{p) has two distinct real roots pi 7 ^ p 2 and po 7 ^ pi, P 2 - Clearly, 
these conditions are eqnivalent to iFi(O) > 0 and po iti The statement (iv) follows 
from fl2.4p . The proofs of the remaining statements follow by direct calcnlations and 
are omitted. □ 

We have the following abont isotropic geodesics. 

Proposition 2.2 Let g G then the following statements hold. 

C. If Pq & at q E then the germ of fl2.ip at (g,Po) is smoothly equivalent to 
p^ = X {Cihrario normal form). 

D. //ord(poll^) = 1 at q E ^ and between the eigenvalues £ 1,62 there are no 
non-trivial resonances 

SiEi + 8262 = Ej, Si,S 2 EZ^, Si + S 2 > 1, j = l,2, (2.5) 

then the germ of cmi at (g,po) is smoothly equivalent to p"^ = y — ex^, where 
e = 6162 7 ^ 0 , A i^Dara-Davydov normal form). The excluded values £ = 0, ^ 
correspond the cases when 6162 = 0 or Si = 82 . 

Z. IfpoW^ at all points on @, then the germ of fl2.ip at (g,po) is smoothly equivalent 
to p^ = y {Clairaut normal form). 

Proof The proof of the case C in the smooth category can be fonnd in [1] (the 
earlier proof in the analytic category can be fonnd in [5]). The proof of the case D 
can be fonnd in [ 6 ] or [7]. The proof of the case Z can be fonnd in [8] (it can be also 
proved nsing the approach in pp for the case C). □ 

Remark 2.2 It is worth observing that Cibrario normal form p^ = x (case C in 
Proposition 12.2p is valid in the analytic category, while Dara-Davydov normal form 
p^ = y — ex'^ (case D in Proposition 12 . 21 ) is not; see [T 8 ] . 
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Table [T] combines the information in Propositions 12.11 and 12.21 and gives a classi¬ 
fication of singularities of the geodesic flow generated by the metric fll.ip . The hrst 
column contains the names of the singularities. The singularities C*, i = 1,2,3 and 
Dg, Dn, Df are topologically stable, while the singularity Z is not and has inhnite codi¬ 
mension. Moreover, the topological equivalence eliminates the modulus parameter e in 
the normal form = y — ex^ in each of three cases Dg, D^, Dj (see mm)- The second 
and third columns are the conditions for the occurrence of the singularities. The last 
column gives the normal forms of the equation of the isotropic geodesics fl 2 .ll) . see also 
Figure [H 


Table 1: Classihcation of singularities of the geodesic flow. 


case 

condition 1 

condition 2 

real roots of M (g, p) 

normal form of (12.11) 

Cl 

Ki{q)<0 

Po ih ^ at g 

Po (simple) 

p^ = X 

C2 

Ki{q)=0 

Po ih ^ at g 

P0i^Pl= P2 

p^ = X 

C3 

Ki{q)>0 

Po iti ^ at g 

Po,pi,P 2 (all simple) 

p^ = a; 

Ds 

Ki{q)>0 

ord (po 1^) = 1 at g 
(12.21) is a saddle at (g,po) on 

Po=Piy^ P 2 

p"‘ = y — ex^ 
e < 0 

Dn 

Ki{q)>0 

ord (po 1^) = 1 at g 
(12.21) has node at {q,po) on 

P0=Ply^ P2 

p"‘ = y — ex'' 
0<e<f^ 

Df 

Ki{q)>0 

ord (po 1^) = 1 at g 
(12.21) is a focus at (g,po) on 

Po = Pi ^ P 2 

P = y — sx 

^ ^ Tfi 

Z 

Ki{q)>0 

Po 1^ at Vg G ^ 

Po = Pi ^ P 2 Vg G ^ 

II 


The following results about the behavior of the geodesics outgoing from a point 
q E ^ with the tangential direction p being a simple root of M{q,p) are established 
in [21]. We start with geodesics passing through a point on the discriminant curve ^ 
with non-isotropic admissible directions. It covers all the cases in Table [D except 6 * 2 • 

Theorem 2.1 f ||21]1 Let Pi, i = 1 or 2, be a non-isotropic simple real root of the 
cubic polynomial M{q,p) at q E S'. Then to the admissible direction pi corresponds a 
unique smooth geodesic passing through the point q. 

The next theorem is proved in [21] for the cases Ci,C^. However, from its proof it 
follows that it is also valid in the case C 2 , since it deals with the root po only. Theorems 
O and O give complete information about the conhguration of geodesics outgoing 
from a point q E S with all possible admissible directions in the cases Ci, Cs. 

Theorem 2.2 f |I21]1 Suppose that the isotropic root po of M{q,p) at q E S is simple. 
Then to pq corresponds a one-parameter family Tq of geodesics outgoing from the point 
q. There exist smooth local coordinates centered at 0 such that the discriminant curve S 
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coincides with the x-axis, the isotropic direction Po{q) = oo and the geodesics 7 ^ G Fq 
are semi-cubic parabolas 

X = ar^Xair), y = r^y„(r), a G (2.6) 

where Xa,Ya are smooth functions, Xq,(0) = 1, Ya{0) = ±1. Here the superscripts 
+ and — distinguish geodesics outgoing from the origin in the semiplanes y > 0 and 
y < t), respectively. The family contains the three types of geodesics: 7 + are 

timelike if a > 1, spacelike if a < 1, isotropic if a = 1; all 7 “ are timelike. See 
Figure |21 




Figure 2: Two examples of configurations of the geodesics outgoing from g G ^ in 
Theorem 12.21 Here the discriminant curve ^ coincides with the x-axis. Timelike, 
spacelike and isotropic geodesics are solid, dashed and bold solid lines, respectively. 


Remark 2.3 1. For any a ^ 0 formula fl2.6p dehnes a semi-cubic parabola, while for 
a = 0 we have the limiting case of the semi-cubic parabolas: the two branches are 
glued together to form the geodesics 7 ^ which are the two semi-axes y > 0 and y < 0. 

2. An interesting problem is to determine whether a geodesic starting from a point 
g G ^ returns to the discriminant curve For instance, two examples of different 
behavior of geodesics are presented in Figure [2J This problem is studied in [22] for 
metrics possessing differentiable groups of symmetries; see also mm for related work. 


3 The main results 

We deal here with the cases Dn, Df,Ds and Z in Table [Hand determine the behavior of 
the geodesics near a point q E S', where the cubic polynomial M{q,p) has the double 
root po = Pi. 















From now on, we will assume that in the cases D^, Dg between the numbers Si, 82 ,1 
there are no non-trivial integer relations 


+ S282 + S3 — £ji 
Sl^l + 52^2 + S 3 = 1 , 


3 

s| := ^ Si >1, Si e Z+. 
1 = 1 


(3.1) 


Clearly, this condition implies the absence of the resonances fl2.5p in Proposition 12.21 
Thus in all cases Df, Dg and Z there exist smooth local coordinates such that 
equation (12.1 p has the normal form = r{x,y), where r{x,y) = y — sx'^. The cases 
Dn,Df,Dg and Z correspond, respectively, to0 <£<^, £>^^^^<0 and £ = 0. 
Then the metric (11. ip has the form 

ds'^ = 0i}{x,y) {r{x,y)dx'^ — dy^) + Q, a;(0,0) = —1, r(x,y) = y — ex^, (3.2) 


where cu is a smooth function and 0 is a smooth metric of the form (II.ip whose 
coefficients are identically zero in the Lorenzian domain. Consequently, its coefficients 
are inhnitely fiat on the discriminant curve r{x,y) = 0. The set of such functions is 
the ideal 911“(r) C (r”) with any n G N (in the ring of smooth functions). 

The assumption ci;(0,0) = —1 can be achieved by multiplication of the metric by 
a non-zero constant, which does not change geodesics. Then y > ex^ (resp. y < ex^) 
is the Lorenzian (resp. Riemannian) domain of the surface. We shall distinguish 
geodesics outgoing from the origin to the region y > ex"^ (resp. y < ex^) using the 
superscripts -|- (resp. —). 

After dividing by appropriate factor (—cu + ...), the held (12.3p is 

2ujrdx + 2pujrdy + Mdp, (3.3) 

with M{x, y,p) = oJxP^ + {ruy + 2uj)p^ — {rux + 2exuj)p — r{rujy -|- cu) - 1 -..., where the 
dots mean terms that belong to the ideal 911“ (r). 

The singular points of the held (13. 3 p are given by the equations r(a:, y) = 0 and 
M{x,y,p) = 0. On the discriminant curve r{x,y) = 0 and away from the origin, the 
cubic polynomial M = p{uJxP^ + 2ijjp — 2exuj) has the simple isotropic root Po = 0 
and two more distinct real roots pi,p 2 being the roots of the quadratic polynomial 
uJxP^ + 2 up — 2exuj whose discriminant is strictly positive on r{x,y) = 0 away from 
the origin. 

Hence in the cases Zl„, Df, Dg at every point on the discriminant curve except the 
origin there are three diherent admissible directions Po,Pi,P 2 , while at the origin two 
of them coincide: po = pi = 0 and p 2 = —2uj/uJx- In the case Z we have Po = Pi = 0 
and P 2 = —2ijjjbJx at all points on the discriminant curve. The case = 0 is also 
included: then M considered as a cubic polynomial over = M U 00 has a simple 
non-isotropic root p 2 = 00 . 

In all cases Dn, Df, Dg and Z, Theorem 12.11 establishes the existence of a unique 
smooth geodesic passing through the origin with non-isotropic admissible direction p 2 . 
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while Theorem 12.21 is not applicable for the double isotropic direction Po = Pi = 0, 
since the spectrum of the linear part of the field fl3.3p at the origin contains three zero 
eigenvalues. 

To study the behavior of the geodesics passing through the origin with the double 
isotropic direction po = Pi = 0 , we consider the blowing up 

T: (x, M,p) (x, i/,p), y = ex^ + up^, u e = M U oo. (3.4) 

The mapping T is one-to-one except on the plane If = {(x,M,p): p = 0}, whose 
image is the line 4/(11) = {(x,i/,p): y — ex^ = p = 0}. The mapping 4/ is a local 
diffeomorphism at all points except on If. It has an inverse defined on \ T(n) given 
by 

T-1 / \ f \ 2 

V {x,y,p}=\x, ,p), r = y-ex. 

Observe that there are no geodesics corresponding to integral curves of the field 
fl3.3p coinciding with the curve 4^(11). Indeed, if 4^(11) is an integral curve of fl3.3p . 
then the identities y = p = 0 hold. Form the first of them we have p = 2ex, which 
contradicts to the second one if e 7 ^ 0. The remaining case e = 0 is more complicated, 
and we give the proof of this statement in Appendix B using the equation of naturally 
parametrized geodesic. 

Away from 4/(11), the map 4/”^ sends the isotropic surface ^ which is an invariant 
surface of the held fl3.3p to the plane u = 1. Hence, the held 03.31) corresponds to a 
smooth held in the (x, M,p)-space (away of H) which has u = 1 as invariant plane. 
Taking this into account, we obtain from 03.3p and 03.4p that the last held, after 
dividing by the common factor p, is 

2ujupdx + (pMi — 2exuj 4-.. .)dp 4- (n — 1) {2uNi + .. .)du, (3.5) 

with 

Mi(x, u,p) = p{upujy 4- Wx)(l — u) + u){2 — u), 

Ni (x, u, p) = up^ujy 4- poJx 4- cn, 

where uj,uJx,(jJy are evaluated at (x,ex^ 4 - np^) and the dots mean terms that belong 
to the ideal 9yt°°{up). 

Since the hrst and second components of 03.5p vanish on the u-axis (x = p = 0), 
the M-axis is the unique integral curve of the held 03.5p passing through any point 
satisfying the conditions x = p = 0 and u{u — l)Ni{x, u,p) 7 ^ 0. From y = ex^ + up^ 
it follows that the map vr o 4/(x,M,p) = (x,p) sends the u-axis in the (x, M,p)-space 
to the origin in the (x, p)-plane. Consequently, the u-axis does not correspond to 
a geodesic. This leads to the necessary condition u{u — 1)A/i(0,m, 0) = 0. Since 
A'i(0, M, 0) = a;(0, 0) 7 ^ 0, we get only two admissible values: u = 0 and u = 1. 

Lemma 3.1 To the admissible value u = 0 does not correspond any geodesic passing 
through the origin in the {x,y)-plane and tangent to the isotropic direction. 
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Proof The germ of the held 03.51) at the origin in the {x,u,p)-spa.ce has the form 
04.7p with n = 3 (Appendix A). Indeed, its components belong to the ideal I generated 
by the components [pMi — 2exu + .. .) and {u — l)(2uNi +.. .), and the center manifold 
W‘^, given by the eqnations pMi — 2exu) = 0 and u = 0, consists of the singnlar points 
of the held 03.5p . The spectrnm of the linear part of 03.51) at any singnlar point is 
(Ai,A 2 , 0 ) with Ai = 2u and A 2 = —2u, so the resonance Ai + A 2 = 0 holds. It 
follows by Theorem 14.41 (Appendix A) that the germ 03.5p is topologically eqnivalent 
to —pdrj. By Theorem 14.71 fAppendix A), it is smoothly orbitally eqnivalent to the 
normal form 04.151) . Moreover, if the 2-jet of 03.51) is generic, it is smoothly orbitally 
eqnivalent to 04.16p . that is, — pdr, + ipdc^. 

Any of these normal forms shows that there are only two integral cnrves passing 
throngh the origin, which coincide with the .^-axis and ? 7 -axis in the normal forms 
coordinates. On the other hand, it is easy to see that the p-axis and w-axis are integral 
curves of the held 03.5p in the initial coordinates {x,u,p). Therefore, the p-axis and 
the M-axis are the only integral curves of the held 03.51) passing through the origin. 
The map vr o T (x, u, p) = (x, y) sends the p-axis and w-axis in the (x, u, p)-space to the 
origin in the {x, p)-plane. □ 

We consider now the admissible variable u = 1 and study the phase portrait of the 
held 03.51) in a neighborhood of its singular point {x,u,p) = (0,1,0). 

Dividing the held 03.5p by the non-vanishing germ {2uNi -|- ...) and making the 
affine change of variable n = m — 1, we transform 03. 5 p into 

pA{x,v,p)dx + {\pBi{x,v,p) - exB2{x,v,p))dp + vdy, (3.6) 

where A and Bi, B 2 are smooth functions such that 

A{x, 0, 0) = Bi(x, 0, 0) = B 2 {x, 0, 0) = 1. 

It is worth observing that 03.6p has an isolated singularity at the origin in the cases 
Ds, Dn, Df (when e ^ 0) and a non-isolated singularity in the case Z (when e = 0). 
We deal with the cases Dg, Dn, Df and Z separately, and start with the case Z. 

3.1 The case Z 

Theorem 3.1 Let S be a surface with a pseudo-Riemannian metric 01 .ip and let 
q ^ S> be a point satisfying the conditions Z in TableUl Then to the isotropic direction 
Po corresponds a one-parameter family Tq of smooth geodesics outgoing from q into the 
Lorenzian domain, while there are no geodesics outgoing from q into the Riemannian 
domain. 

There exist smooth local coordinates centered at the origin such that the metric has 
the form 03.2p and the geodesics 7 + G Tq outgoing from the origin are 

y =-^x'^ + ax'^il+ Ya{x)), 1^(0) = 0, a e M. (3.7) 
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The family 03.71) contains all three types of geodesics: timelike if a > t], spacelike if 
a < 0 and isotropic if a = 0 ; see Figure [21 right. 


Proof The germ 03.61) with e = 0 at the origin in the {x, v, p)-space has the form 04.71) 
with n = 3 (Appendix A). Indeed, its components belong to the ideal / = {pA,v) = 
{p,v) and the center manifold W‘^, given by the equations p = v = 0, consists of its 
singular points. At any singular point sufficiently close to the origin, the spectrum of 
the linear part of 03.6p is (1,|,0). Then, by Theorem 14.51 (Appendix A), the germ 
03. 6 p at the origin is orbitally smoothly equivalent to 

(2^ + T{C,W)di + Vdr,- (3.8) 

The restriction of the field 03.8p to each invariant leaf ( = const is a resonant node 
with eigenvalues 2 : 1 and associated eigenvectors <9^, 2dx + dp. It has an inhnite family 
$ of integral curves all with the common tangent direction 2dx + dp and a unique 
integral curve with the tangent direction dy, which coincides with the n-axis and does 
not correspond to a geodesic. 

By Lemma [4.21 (Appendix A), to prove that (p{C) vanish on each leaf ( = const, it 
is sufficient to produce a smooth integral curve of the family <h. We use the integral 
curve X = 2p, V = 0 corresponding to the isotropic geodesic y = which is a solution 
of equation p^ = y. Hence (^(C) = 0, and the normal form fl3.8p becomes 

2^d^ + T]dp. (3.9) 

Comparing fl3.6p with e = 0 and 03.81) . one can show that the conjugating diffeo- 
morphism {x,v,p) {^,vX) can be chosen in the form 

^ = T] = pf{x,v,p), C = X - 2p + g{x,v,p), /(O) = 1, geTl^, 

where 971^, /c > 0 , is the ideal of /c-flat functions in the ring of smooth functions. 

The field 03.9p has the invariant foliation = cr]"^, c G = M U oo. Returning 
to the coordinates {x,v,p), we get the invariant foliation of the held 03.6p given by 

V = cp^f{x,v,p), /(O) = 1. (3.10) 

The value c = oo corresponds to the invariant plane H = {p = 0} which does not 
give any geodesics. For c G M, solving O3.10p in n in a neighborhood of the origin gives 
V = cp‘^'ipc{x,p), where 'ifcix,p) is a smooth function with '0c(O,O) = 1. Returning to 
the initial variables {x, y,p) by mean of the map T (formula 03.41) 1. we get the invariant 
foliation of the held 03.3p : 

y = p'^{l + cp‘^if^{x,p)), '0c(O,O) = l, (3.11) 

see Figure 12] (left). The value c = 0 corresponds to the isotropic surface = y. 
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Observe that all invariant leaves of the foliation fl3.1ip belong to the semiplane 
y > 0. The plane y = 0, which does not belong to the foliation fid.lip . is an invariant 
plane of the held fl3.3p foliated by the vertical lines, which are integral curves of fl3.3p . 

After dividing by the common factor the restriction of the held fl3.6p to each 
invariant leaf fl3.1ip is given by 


dp uj — cp'^'ipcipJ + P^x + y^y) 
dx 2u}(l + cp'^'ipc) 


1 

2 


cp^(l + hc{x,p)), 


(3,12) 


where w, oOy are the restrictions of cn, ooy to the given leaf. In a neighborhood 
of the origin, fl3.12p dehnes a smooth and regular direction held, and every leaf 03.111) 
contains a unique smooth integral curve of the held 03.6p passing through the origin 
in the (x, |/,p)-space with the tangent direction dx, see Figure [3] (center). Projecting 
these integral curves to the (x,|/)-plane (equivalently, integrating equation 03.12p 'l. we 
get a one-parameter family of geodesics in the statement of the theorem with a = 
see Figure [3] (right). □ 



A 


y 


X 

> 



Figure 3: The case Z: invariant foliation 03.lip of the held 03. 3 p (left), integral curves of 
03. 3 p on an invariant leaf (center) and geodesics 7 ^ G Fq (right). Here the discriminant 
curve ^ coincides with the x-axis. Timelike, spacelike and isotropic geodesics are solid, 
dashed and bold solid lines, respectively. 

To illustrate Theorem 13.11 we consider the simplest example ds^ = dy^ — ydx^. In 
this case, equation 02.31) can be studied using qualitative methods, see, for example, [22] 
(Section 3). The Lagrangian of the length functional L = \/p'^ — y does not depend 
on the variable x, hence the held 02.3p possesses the energy integral H = L — pLp. 
After evident transformations, equation H = const can be reduced to 

p^ = y — ay‘^, a G M, (3.13) 

which is a family of implicit diherential equations of Clairaut type. 
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Every (unparametrized) geodesic in the metric = dy^ — ydx^ is a solution of 
equation fl3.13p . Conversely, every solution of fl3.13p is a geodesic except the horizontal 
lines y = const, each of which is the envelop of the family of integral curves of fl3.13p 
for a given a (see IZ21), For instance, the value a = 0 corresponds to the isotropic 
surface = y {a. parabolic cylinder) and gives the isotropic geodesic y = \x^ passing 
through the origin. For determining non-isotropic geodesics, observe that every invari¬ 
ant surface fl3.13p is a cylinder whose generatrices are parallel to the x-axis and the 
base is an ellipse (if a > 0) or a hyperbola (if a < 0). In the last case, the hyperbolic 
cylinder = y — ay"^ consists of two connected components: positive and negative 
lying in the domains |/ > 0 and y < a~^, respectively. See Figure IH left. 

Positive components of the hyperbolic cylinders {a < 0) together with all other 
cylinders {a > 0) form an invariant foliation of the form 03.101) (Figure IU left). Neg¬ 
ative components of the hyperbolic cylinders do not intersect the plane y = 0, and 
consequently, they do not contain integral curves whose projections to the (x, y)- 
plane are geodesics passing through the x-axis. Moreover, from 03.13p it follows that 
u = y/p"^ = 1/(1 — ay), which shows that u —?■ 1 along any geodesic tending to the x- 
axis with isotropic tangential direction. This gives a graphical explanation why u = 1 
is the only admissible value of the variable u ai y = p = t). 

Thus to every a > 0 corresponds a geodesic 7^ G Tq which is timelike if a > 0 or 
isotropic if a = 0. To every a < 0 corresponds a spacelike geodesic 7 + G Tq, whose lift 
belongs to the positive component of the hyperbolic cylinder p^ = y — ay"^. In contrast 
to this, the negative component of the same cylinder is hlled with integral curves of 
the held fl2.3p whose projections on the (x, |/)-plane are separated from the x-axis by 
the horizontal strip < y <0. See Figure Hj right. 




Figure 4: The invariant foliation p^ = y — ay'^ in the (x, ?/,p)-space (left) and the 
corresponding geodesics (right). Timelike, spacelike and isotropic geodesics are solid, 
dashed and bold solid lines, respectively. 
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3.2 The cases Ds,Dn,Df 

We take the metric as in fl3.2p with e 0,^ near the point q E The germ fl3.6p 
with e: 7 ^ 0 at the origin is hyperbolic with the spectrum ( 1 , ei, £ 2 ), where (ei, £ 2 ) is the 
spectrum of the field fl2.2p at the origin, so £i +£2 = | and 8162 = e. The corresponding 
eigenvectors are d^, + Sidp and + £ 2 ^p- 

Recall that we assumed that, in the cases Dn,Ds, there are no non-trivial reso¬ 
nances fl3.ip between £ 1 , 62 ,1- Hence, by Theorem 14.11 (Appendix A), the germ 03.61) 
at the origin is smoothly equivalent to 


+ eir]dp + e 2 Cdi;. (3.14) 

Remark 3.1 It follows by comparing 03. 6 p and 03.141) that the conjugating diffeomor- 
phism (x,v,p) I—)■ and its inverse can be chosen in the form 

X = 77(1-f (pi) + C (1 + 952), p = ? 7 (£i(ps) + C(£2 + 754), V = ^; 

_ x{e 2 + iJi) - p{l + i> 2 ) x{ei + iJs) - p{l + ^> 4 ) ^ (3-15) 

h ) s ) s n, 

82 — £1 £1 — £2 

where p* = p, () and ipi = 'ipi{x, v,p) are smooth functions vanishing at the origin. 


Theorem 3.2 Let S be a surface with a pseudo-Riemannian metric HHI and let 
q E 3> be a point satisfying the conditions Dg in Table [3 Then to the isotropic 
direction pq corresponds a one-parameter family Tq of C'^-smooth geodesics outgoing 
from q into the Lorenzian domain, while there are no geodesics outgoing from q into 
the Riemannian domain. 

There exist smooth local coordinates centered at the origin such that the metric has 
the form 03.2p and the geodesics 7 + G Tq outgoing from the origin are 

y = Yaix) = o(x^), « G M, (3.16) 

together with the isotropic geodesic 

p = y(x), y(x) = o(x^), £2 < 0. (3.17) 

The family 03.16p contains all three types of geodesics: timelike if a < Q, spacelike if 
a > 0 and isotropic if a = 0; see Figured right. 

Proof By Lemma ITTI fAppendix A), equation + C 2 |CI^^^^ Tcs'C = 0 with any 

constants q defines an invariant surface of the held 03.141) . Since £ < 0, the eigenvalues 
£i and £2 have different signs. We assume, without loss of generality, that £i> \ and 
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82 < 0. Hence, the invariant surface Ci\r]\^^^^ + C 2 |C|^^^^ + = 0 passes through the 

origin if and only if Ci = C 3 = 0 or C 2 = 0. This yields the invariant foliation 

^ = a\r]\^^^\ aeMP^ = MUoo, (3.18) 

where a = 00 corresponds to the plane 77 = 0 being the limit set of the cylindric 
surfaces ^ as a —)■ cxo. The restriction of the held fl3.14l) to each invariant 

leaf 03.181) is a saddle with the separatrices 

{^ = "^ = 0} and C = 0} for « 7 ^ 00 , 

{■C = "^ = 0} and {77 = C = 0} for a = 00 , 

where the separatrix {.^ = 77 = 0} is the intersection of all invariant leaves 03.18p . The 
corresponding invariant foliation of the held 03.3p is as in Figure O (left). 

Taking into account 03.151) . u = v + 1 and y = ex^ + one can see that the 
separatrix { 7 / = C = 0} corresponds to the n-axis in the (x, 7;,p)-space and to the 
origin in the (x,7/,p)-space. The remaining separatrices = «|?7|^/^hC = 0} 

{^ = h = 0} yield the family of integral curves of the held 03.51) passing through the 
point ( 0 , 1 , 0 ) in the (x, M,p)-space: 

{x = 77 (l+(pi,„( 77 )), p = 77(61 + ^sAv)), u = l + a e R, (3.19) 

{x = C(1 + ^ 2 ( 0 ), P = Ci ^2 + AiO), u = 1 }, (3.20) 

where AAv) = PiA\vA‘'\vA) for 7 = 1,3, and AiO = V5i(0,0,C) for 7 = 2,4. 
Observe that the integral curve 03.19p with a = 0 and the integral curve 03.201) belong 
to the isotropic surface {77 = 1}. 

Since the right-hand sides of x = 7/(1 -1- ^i,a( 77 )) and x = C(1 + ^ 2 ( 0 ) ^ire and 
0°°-smooth respectively, both equations can be solved for 77 . Hence, from 03.191) and 
O3.20p one can express the variables p and u as functions of x. To complete the proof, 
observe that the map tt o v]>(x, u,p) = (x, y) sends the family of integral curves 03.191) 
to the family of geodesics 03.161) and the integral curve 03.201) to the geodesic 03.17P : 
see Figure E] (right). □ 

Remark 3.2 In the case Z the family Fq given by 03.71) contains a unique isotropic 
geodesic y = |x^ (the bold line in Figure El right) which separates Fq into timelike 
and spacelike geodesics. In the case Dg the family Fq contains two isotropic geodesics, 
corresponding to the integral curve 03.19p with a = 0 and the integral curve O3.20p . 
The hrst (the bold line in Figure [5l right) separates Fq into timelike and spacelike 
geodesics, while the second (the double line in Figure^ right) does not. The difference 
between the families Fq in the cases Z and Dg is related to topological structure of the 
invariant foliations of the held fl3.3p in the cases Z and Vg; compare Figures [3] (left) 
and [5] (left). In the case Z all invariant leaves intersect on the line 4/(n) only, while 
in the case Vg they intersect on the dotted line 4^(n) and on the double line which 
corresponds to {.^ = 77 = 0} in the normal form fl3.14p . 
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Figure 5: The case Dg'. invariant foliation of the held fl3.3p (left), integral curves of 
fl3.3p on an invariant leaf (center) and geodesics 7 ^ G Fq (right). On the right: timelike 
and spacelike geodesics are solid and dashed lines, respectively. The bold line and the 
double line are isotropic geodesics, the dotted line is the discriminant curve. 


Theorem 3.3 Let S be a surface with a pseudo-Riemannian metric uni and let 
q E ^ be a point satisfying the conditions or Df in TableUl Then to the isotropic 
direction pq corresponds a two-parameter family Fq of C"^-smooth geodesics outgoing 
from q into the Lorenzian domain, while there are no geodesics outgoing from q into 
the Riemannian domain. 

There exist smooth local coordinates centered at the origin such that the metric has 
the form fl3.2p and the field fl3.3p has the invariant foliation 

y = ex^p'^{laYa{x,p)), Ya a G M, (3-21) 

see Figure 0 left. The restriction of the field fl3.3p to every invariant leaf (13.2 ip is 
a node {for Dn) or a focus {for Df) with spectrum 61 , 62 , and the geodesics 7 ^^ G Fq 
are projections of the integral curves of (13.3p on the leaves (I3.2ip to the {x,y)-plane 
(Figure [HI right). The geodesics 7 ^^ G Fq are timelike if a <t], spacelike if a > and 
isotropic if a = 0 . 


Proof We start with the question of linearizability of the germ (13. 6 p at the origin. 
In the case D^, the germ (13.6p is smoothly equivalent to the linear normal form (I3.14p 
with 0 < £ 1 , £2 < I (Theorem [TT] in Appendix A), and the conjugating diffeomorphism 
has the form (13.151) . 

In the case Df, Theorem 14.11 does not apply because the spectrum (l,£i,£ 2 ) with 

61 = a + bi, 62 = a — ib, a = |, b = i\/l 6 £ — 1, 6 = 6162 > (3.22) 

has the resonance 2(£i + £ 2 ) = 1 of the order |s| = 4 impeding the linearizability in the 
CAsmooth category. Consider along with (13.ip the resonances obtained by taking the 
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Figure 6: The cases Dn and Df-. invariant foliation of the field fl3.3p (left), integral 
curves of the held fl3.3p on an invariant leaf and their projection to the (x, |/)-plane in 
the cases (center) and Df (right). 


real part of both sides of fl3.ip . The spectrum fl3.22p has hve resonances of this type 
of order greater than one: Re(siei + 8262 ) = 1, |s| = Si + S 2 = 4. Thus in the case Df, 
the germ fl3.6p satishes the condition of Theorem 14.31 (Appendix A) with k = 1, and 
consequently, it is C^-smoothly equivalent to 

+ {ap + bQdr, + (aC - for/)«9^ (3.23) 


with a, b as in fl3.22p . 

Thus in the case Dn (resp. Df) the germ fl3.6p is C°° (resp. C^) smoothly equivalent 
to the linear normal form fl3.14p (resp. fl3.23p h and the conjugating diffeomorphism 
has the form 

x = ri^i + CP2, P = VP3 + CP4, v = 
r] = xV'i + pV' 2 , C = x'ljjs + pV^4, ^ = v, 

for some C°° (resp. C^) smooth functions tpi = and X = '^i(x, n,p). 

By Lemma l4T] (Appendix A), equations 

f = a(|r,r'''' + icr-'”), (3.25) 

5 = a(y + e)\ (3.26) 

where a G MP^, define invariant foliations of the linear helds fl3.14p and fl3.23p . re¬ 
spectively. Here a = 00 gives the .^-axis in the normal coordinates, i.e., the n-axis in 
the (x, n,p)-space, and consequently, the origin in the (x, |/,p)-space. Hence we shall 
consider both equations fl3.25p and (13.261) with a G M only. 

Substituting expressions fl3.24p for ^,r],C in (13.251) and (I3.26p . gives the invariant 
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foliation of the field fl3.6p in the case Dn and Df respectively: 

V = + (3.27) 

V = a{\x'ijji +P'iIj2\‘^ + \x^p3 +p-ipi\'^Y. (3.28) 

Note that the right-hand sides of fl3.27p and fl3.28l) are at least C^-smooth (in the 
case Dn, 0 < £ 1,62 < |)- By the implicit function theorem, fl3.27p and fl3.28p can 
be solved in v near the origin. In both cases we get the equation v = aYa(x,p) with 
a C^-smooth function vanishing at the origin. Substituting this expression into 
y = ex'^ + (1 + gives the invariant foliation fl3.2ip of the held fl3.3p . 

Similarly to the case Z, one can divide the restriction of the held fl3.3p to every 
invariant leaf by the common factor p (see formula fl3.12p h Then the restriction of the 
held fl3.3p to every invariant leaf fl3.2ip becomes a node (in the case Dn) and a focus 
(in the case Df) with the spectrum { 81 , 62 ) (Figure [ 6 l center, right). □ 

Figure [7] presents computer generated pictures of geodesics on a surface endowed 
with the metric fl3.2p with a; = — 1 and 0 = 0 in the cases Dg, Dn, Df. 



Figure 7: Geodesics (solid lines) in the metric fl3.2p with uj = —1 and 0 = 0 in the 
cases Dg (left), Dn (center), Df (right). The dotted line is the discriminant curve. 


4 Appendix A. Local normal forms of vector fields 

Here we give a brief survey of local normal forms for vector helds in real phase space, 
which were used in this paper (for more details, see also surveys in [21 [3]). All vector 
helds are supposed to be smooth (that means unless stated otherwise). By 
we denote the class of analytic mappings. For convenience, we also present vector 
helds as autonomous diherential equations, where the diherentiation is by an auxiliary 
parameter playing a role of time. 
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Definition 4.1 Two vector fields Vi and V 2 are C'^-smoothly (resp. topologically) 
equivalent, if there exists a C'^-diffeomorphism (resp. homeomorphisni) h: M” 

that conjugates their phase flows and g^, i.e., h o g\ = g\ o h. Here k is an integer 
number (finite-smooth equivalence) or 00 (infinite-smooth equivalence) or uj (analytic 
equivalence). 

Definition 4.2 Two vector fields Vi and V 2 are orbitally C^-smoothly (resp. topo¬ 
logically) equivalent, if there exists a C'^-diffeomorphism (resp. homeomorphism) that 
conjugates their integral curves (orbits of their phase flows). 

Remark 4.1 The second definition slightly differs from the generally accepted defini¬ 
tion of the orbital equivalence, where coincidence of the orientation of integral curves 
is also required. In fact, our definition is naturally related to directions fields, whose 
integral curves do not have a orientation a priori. 

A great deal of work is done to bring a vector field to a normal form in a neigh¬ 
borhood of its singular point under a chosen equivalence relation. In particular, the 
germ of a vector field at its singular point is called linearizable (in a certain category) 
if it is equivalent (in this category) to its linear part. Of course, not all vector fields 
are linearizable, and when it is not, the question is what normal formal can we get. 

Definition 4.3 Let V be the germ of a vector field with singular point the origin 
0 in M"" and A = (Ai,...,A„) E C” be the spectrum of V at 0. The germ V is 
called hyperbolic if the spectrum A does not contain neither zeros nor pure imaginary 
eigenvalues. The germ V is called partially hyperbolic if the spectrum A contains at 
least one eigenvalue with non-zero real part. 

4.1 Hyperbolic germs 

Let V be the germ of a hyperbolic vector field at 0. For topological equivalence the 
only local invariants are the scalars = sgn(Re Aj). The Grobman-Hartman Theorem 
Enn] states that V is topologically equivalent to the field = 0 *^*, i = 1,... ,n. 

For the G^-smooth equivalence, resonances of two types play an important role: 

Aj-(s, A) = 0, SiEZ+, je {!,..., n}, (4.1) 

Re(Aj - (s. A)) = 0, Si e Z+, j E {1,.. . ,n}, (4.2) 

where (s. A) = SiAi s„A„ is the standard scalar product. In both relations fl4.ip . 

fl4.2p the natural number |s| = Si-|- - • - -l-Sn is called the order oi the resonance, and = 

called the resonant monomial. In general, non-trivial resonances fl4.ip . 
fl4.2p of orders |s| > 2 are obstacles for the germ V to be G^-smoothly linearizable. 
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Theorem 4.1 ([21 SI 129^ If the spectrum A does not have resonances fl4.ip of any 
order |s| > 2, the germ V is C°°-smoothly linearizable. Moreover, z/A G M" and has 
only trivial resonances (BUI, the germ V is C°°-smoothly eguivalent to 


(4.3) 


* 1, . . . , 77.. 


Theorem 4.2 ([21 SI E]) If \ E M”, then the germ V is C°^-smoothly eguivalent to 
the field 



(4.4) 


where ats ^ 0 only if the resonance fl4.ip holds. 

In other words, Theorems 14.11 and 14.21 states that the normal forms fl4.3p . fl4.4p 
contain resonant monomials only. Here the linear terms Xifi correspond to trivial 
resonances (03, while all remaining terms correspond to non-trivial resonances 
(14.1 p if they exist. Under some additional restrictions, Theorems 14.11 14721 are valid in 
C‘^ category. For instance, in the case A G M” it is sufficient to require that all Aj have 
the same sign. However, if Aj have different signs, the restrictions are much stronger, 
especially if A has non-trivial resonances fTheorem 14.2p . see [21 E] and the references 
therein. 

On the contrary, requirements become weaker if we consider U^-smooth equivalence 
with k < oo, see e.g. [21 |25l [26]. In this paper, we need the following result. 

Theorem 4.3 ([21 127] 1 Suppose that for some positive integer k the spectrum A has 
neither resonances (03 of order 2 < |s| <2k nor resonances (14.2p of order |s| = 2k. 
Then V is -smoothly linearizable. 

We establish below the existence of invariant foliations of a hyperbolic linear vector 
held in 3-dimensional real phase. Assume, without loss of generality, that one of the 
eigenvalues Aj of is equal to 1. It is sufficient to consider the following two helds 


— Ai^i, ^2 — A26, ^3 — 6) 


6 = «i+/56), 6 = (<2-/56), 6 = 6, 


(4.5) 

(4.6) 


where Ai, A 2 and a, (3 are real and non-zero. 

Lemma 4.1 For any real constants Ci the eguations 


cii6r/^^+c2i6r/^=+c36 = o, 
ci(e? + e2)'/'“ + c36 = o. 


define invariant surfaces of the fields (14. 5 p and (14.61b respectively, 
The proof is trivial and is omitted. 
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4.2 Partially hyperbolic germs 

Let W^, W^, be the unstable, stable and center manifold of the partial hyperbolic 
germ V at 0, and let di = dimbh*, i = s,u,c. Set d = ds + du, then d,dc > 0 and 
d + dc = n. One can choose local coordinates (^i,.. .,G ■ ■ ■ ,^d) £ 

C = (Ci,...,a)efLL 

Theorem 4.4 ([21 [12]) The germ V is topologically eguivalent to the direct product 
of d-dimensional standard saddle (the first d eguations) and the restriction ofV to the 
center manifold (the last dc eguations): 

i !)•••) dg, i dg -|- 1,..., d, 

(j = Zj(fifi j = 1,..., dc- 


In this paper, we deal with a special class of partially hyperbolic vector helds, 
which were studied by many authors, see e.g. [231I3I]. From now on, we assume that 
all components of the germ V belong to the ideal I (in the ring of smooth functions 
vanishing at 0 ) generated by two of them. 

More specihcally, such a germ V has the form 

i = v, T] = w, Cj = ajV + l3jW, j = l,...,n-2, (4.7) 

where v,w and are smooth functions of the variables (,,ri,(j- The components 

of the germ 04.711 belong to the ideal / = {v,w), and the spectrum of V contains at 
most two non-zero eigenvalues: A = (Ai, A 2 , 0,..., 0). 

We shall further assume that Re Ai 7 ^ 0 and Re A 2 0. Hence the center manifold 
IF^ = {u = tn = 0} C is a smooth manifold of codimension 2 and the restriction 
of the held V to is identically zero, so consists of singular points of V. By 
Theorem 14.41 the germ V is topologically equivalent to 

i = V = a2V, 6 = 0) J = 1, • • •, ^ - 2, 

where = sgn(Re \i), i = 1, 2. 

For C^-smooth classihcation of the germ fl4.7p . we need to introduce two types 
of resonances between the non-zero eigenvalues Ai,A 2 being, in fact, partial cases of 

(El): 


siAi -|- S 2 A 2 — 0, Si G Z_|_, i — 1,2, (4-8) 

SiAi-|-S 2 A 2 = Aj, Sj G Z_|_, i,j = 1,2. (4-9) 

For simplifying the presentation, further we shall always assume that |Ai| > IA 2 I 
and exclude from consideration trivial resonances fl4.8p (si = S 2 = 0) and fl4.9p (si = 1, 
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j = 1 or S 2 = 1, j = 2). Then the absence of resonances fl4.9p implies the absence of 
fl4.8p . On the other hand, in the absence of 04.81) . resonances 04.91) may have only the 
form Ai = mX 2 , for positive integers m. 

Given the germ 04.7p with Re Ai 7 ^ 0 and Re A 2 7 ^ 0 we choose local coordinates 
^, 1 ] (called hyperbolic variables) and ( = {(i,... Xn- 2 ) ^ (called non-hyperbolic 
variables) snch that the ideal / = X,ri), and conseqnently, the center manifold is 
given by ^ = rj = 0. The linearization of V with respect to the hyperbolic variables 
has two eigenvalues which are continuous functions Ai(C) and A 2 (C) of C G W^. We 
have Aj(0) = Xj, j = 1,2 and Xj as above. 

An analogue of Theorems 14.11 and 14.21 is the following. 


Theorem 4.5 (HOI [23]) Suppose that between Ai(C) and A 2 (C) there are no non¬ 


trivial resonances of any order |s| > 1 for all ( sufficiently close to zero. Then 
the germ (I4TI) is C°°-smoothly equivalent to 

i = X, r] = Y, Cj = J = 1, • • •, - 2, (4.10) 

where X,Y are smooth functions of such that the ideal I = {X,Y) = {f,ri). 

Moreover, if in addition the eigenvalues Ai,A 2 G M, then 

X = XrX)i + XC)r}^, Y = X2 {C)v, (4.11) 

where v?(C) ^ 0 only if Ai = mX 2 with some natural m>l. 


Remark 4.2 If the pair (Ai,A 2 ) belongs to the Poincare domain (i.e., Ai and A 2 are 
real and of the same sign or complex conjugate), then the condition 

5iAi(C) + 52A2(C)^0, gZ+, * = 1,2, VCgW^ (4.12) 

follows from 04.81) . Moreover, in this case Theorem 14.51 is valid in category. How¬ 
ever, if the pair (Ai,A 2 ) belongs to the Siegel domain (i.e., Ai and A 2 are real and 
of different signs), the condition 04.12p is equivalent to Ai(C) : A 2 (C) = const for all 
C G WX 


The conditions in Theorem 14.51 become weaker if we consider G^-smooth equiva¬ 
lence with k < 00 . Set 


N{k) = 2 


{2k+ 1) 


max I Re Ai ^2 


min I Re Ai ^2 

where the square brackets is the integer part of a number. 


+ 2, keN, 


Theorem 4.6 f [ jlOp For any k E N, the statements in Theorem still hold 
true if C°° is replaced with and the inequalities 1 < |s|, 1 < m are replaced with 
1 < |s| < N{k), 1 < m < N{k) respectively. 
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The normal form fl4.10p . fl4.1ip can be fnrther simplified. For our purposes, we are 
interested in orbital normal form in the case when the resonance Ai(C) = 'm^ 2 {C) holds 
at all C G W^. Then, dividing by A 2 (C); from fl4.10p . fl4.1ip we get the orbital normal 
form 

i = + V = V, 0 = 0) j = l,...,n-2. (4.13) 

where the smooth functions (p{C) and ' 0 (C) vanish simultaneously. 

The following lemma gives a simple geometric criterion for ' 0 (C) = 0 , which is 
important for applications. 


Lemma 4.2 Let V be the germ of a field from Theorem \4-5\ with the normal form 
fl4.10p . 04.1111 and the resonance Ai(C) = 'mA 2 (C); m > 1, holds at all points ( G W^. 
Then in the orbital normal form 04.13p . 0(C) = 0 if and only if V has a C^-smooth 
integral curve that passes through the corresponding point ( G with the tangential 
direction parallel to the eigenvector with A 2 (C)- 


Proof The held 04.13p can be integrated explicitly. It has the invariant foliation 
C = const and each leaf contains a single integral curve rj = 0 with tangential direction 
and one-parameter family of integral curves 

f = ri"^{c + fi>{()\n\ri\), c = const, (4-14) 

with the common tangential direction 0^ at the point f = rj = 0. All the curves 04.1411 
are C^^^-smooth at ^ = -z] = 0 if (^(C) 0 0 and C'°°-smooth at zero if (^(C) = 0 . 

Given ( G W^, the existence of at least one G^-smooth integral curve passing 
through the point ^ = 'z; = 0 (the intersection of with the corresponding invariant 
leaf C = const) with the tangential direction parallel to the eigenvector with A 2 (C) is 
equivalent to the condition 0 ( 0 ) = 0 . □ 

Theorem 4.7 f [lOp I23jl Suppose that the resonance Ai(C) + A 2 (C) = 0 holds at all 
singular points C G and ReAi 0 0, Re A 2 0 0. Then, for any natural k, the germ 
V is -smoothly eguivalent to 

C = C(Ai(C)+ p4>i(p,C)), T] = r]{X2iC) + P^2ipX)), 

(4.15) 

0 =pT 0 p,C), j = l,...,n- 2 , 

where 4’i(p, C) and T0p, Q are polynomials in p = fp of degrees N{k) — 1. 

If 4^00, 0) 0 0 for at least one j = 1,... ,n — 2, then the germ V is C°°-smoothly 
orbitally eguivalent to 

C = C, V = -V, C = C7, j = l,...,n- 2 . (4.16) 


Theorem 14.71 is not valid in C‘^ category. 
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5 Appendix B. Naturally parametrized geodesics 

Naturally parametrized geodesics can be defined as extremals of the action functional 

J{l) = j + cf) dt, y = 

7 

where 7 C S' is a differentiable curve. The corresponding Euler-Lagrange equation 
reads 

I 2{ax + by) = (c^ - 2by)y‘^ - 2ayXy - a^x^, 

I 2{bx + cy) = {ay - 2b^)x^ - 2c^xy - Cyij^. 

The definition of geodesics as auto-parallel curves in the Levi-Civita connection gen¬ 
erated by the metric fll.ip leads to the same equation fIS.ip . Equation fIS.ip defines a 
direction field on the tangent bundle TS. The standard projectivization TS ^ PTS 
sends this direction field to the field parallel to fl2.3p . see [ 22 ] , 

Firstly, using equation fIS.ip of parametrized geodesics, we prove the omitted state¬ 
ment in the case Z that the line T(n) = {y = p = 0} does not correspond to a 
geodesic. Recall that in the case Z there exist local coordinates such that 

ds^ = {yu + .. .)dx^ -|- (0 . .)dxdy — (cn -|- .. djdy"^, a;( 0 , 0 ) = — 1 . 

where the dots mean terms that belong to the ideal DJl°^{y). 

Using appropriate change of variables y i-7 yu{x,y), where n is a solution of equa¬ 
tion cUx + b/y = 0 with the condition n( 0 , 0 ) 7 ^ 0 , one can bring locally the metric 
to the diagonal form ds^ = adx'^ -|- cdy"^ with the coefficients a = yuu + ... and 
c = —uj{u -|- yUy). Substituting y = b = 0 into the second equation of fIS.ip . we get 
ay{x, 0)i;^ = 0. Since 0 ^( 0 , 0) 7 ^ 0, this yields x = 0, and the restriction of the system 
(EHi to p = 0 has only constant solutions, which are not geodesics. 

At first sight it contradicts the fact established in [10]: ^ is an invariant surface 
of the field fl2.3p . and consequently, any trajectory of fl2.3p that lies entirely in ^ 
after the projection on the (x,|/)-plane gives a geodesic or a point. (Example: for the 
metric ydx"^ — dy^ the isotropic surface = p is filled with one-parameter family of 
integral curves intersecting T(n) transversally. Projecting this family down, we get 
the isotropic geodesics y = ^{x — c)^.) In fact, there is no contradiction: the curve 
\h(n) C ^ consists of singular points of the field fl2.3p . and every such a point is a 
trajectory of fl2.3l) . 

Consider the family Tq of geodesics outgoing from a point g G ^ with the isotropic 
direction po. Choose the natural parametrization so that the motion along geodesics 
proceeds toward q. In the paper [22], it was proved that in the cases Ci,C '3 any 
geodesic 7 G Tq incomes into the point q in finite time with infinite velocity. The same 
statement is valid in the case C 2 , since it deals with the root po only. 
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In the case Z, the same result follows from the asymptotic formula established 
in Theorem 15.11 below. The cases Ds,D„, can be considered similarly, the case Df is 
excluded from consideration, since every geodesic 7 G Tq intersects the discriminant 
curve S' inhnite number of times in any neighborhood of q. 

Theorem 5.1 The natural parametrization of geodesics dSTD is given by the formula 
a; = fs (1 + Xait^)), where Xa{-) are smooth functions vanishing at zero. 

Proof Choosing the local coordinates in Theorem 13.11 from the formula 03.711 we 
have y = + 0 (a:^), ij = {\x + 0{x^))x, and y = {\x + 0{x‘^))x + (| + 0{x‘^))xf‘. 

Substituting these expressions together with the coefficients a, b, c from 03.2p into the 
hrst equation in 05.ip . after a straightforward transformation we obtain 

^ = (--+ fa{x))x, (5.2) 

where faix) are smooth functions. 

Equation 05.2p dehnes the natural parametrization uniquely up to non-degenerate 
affine transformations of the f-axis. Integrating it, we get In |x| = —2 In |a:| +Fa{x) + C, 
and X = where Fa is the primitive of fa- Without loss of generality put 

-^o(O) = 0 and K = ^ (this corresponds to the choice of the initial velocity of motion 
along the geodesic). Then we arrive at the differential equation ^ whose 

general solution is f = a:^(l -|- Ta{x)) -|- to, where Ta{x) is a smooth function vanishing 
at zero. Setting to = 0 and inverting, we get x = fs (l -|- Xait^)). □ 

As an example, for the metric = dy"^ — ydx"^ the system 05.1 p reads yx = —xy, 
2y = —if'. Substituting here the isotropic geodesic y = we get x = k{t — to) 3 . 
Substituting ?/ = 0 (the line T(n)), we get i = 0, that is, y = 0 is not a geodesic as 
was stated above. 
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